Nonequilibrium quench dynamics in quantum 
quasicrystals 



Ferenc Igloi^^^, Gergo Roosz^'^ 

^ Wigncr Research Centre, Institute for Solid State Physics and Optics, H-1525 
Budapest, P.O.Box 49, Hungary 

^ Institute of Theoretical Physics, Szeged University, H-6720 Szeged, Hungary 

E-mail: igloi . ferencSwigner.mta.hu 

E-mail: gergo_roosz@t itan . physx . u-szeged . hu 

Yu- Cheng Lin 

Graduate Institute of Applied Physics, National Chengchi University, Taipei, 
Taiwan 

E-mail: yc . linOnccu. edu. tw 

Abstract. We study the nonequilibrium dynamics of a quasiperiodic quantum 
Ising chain after a sudden change in the strength of the transverse field at zero 
temperature. In particular we consider the dynamics of the entanglement entropy 
and the relaxation of the magnetization. The entanglement entropy increases 
with time as a power-law, and the magnetization is found to exhibit stretched- 
exponential relaxation. These behaviors are explained in terms of anomalously 
diffusing quasiparticles, which are studied in a wave packet approach. The 
nonequilibrium magnetization is shown to have a dynamical phase transition. 
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1. Introduction 

Recent experimental progress in ultracold atomic gases in optical lattices [H El El 
El El m [7l El El HD] has opened fascinating new perspectives on research in the 
field of isolated quantum systems, both in equilibrium and out of equilibrium. In 
experiments the form of atomic interactions can be suddenly changed by tuning an 
applied magnetic field near a Feshbach resonance, which is known as a global quantum 
quench. On the theoretical side, one is interested in the time-evolution of different 
observables, such as the order parameter or some correlation function, after a quench. 
Fundamental questions concerning quantum quenches include (i) the functional form 
of the relaxation process in early times, and (ii) the properties of the stationary state 
of the system after a sufficiently long time. 

Many results for quantum quenches have been obtained for homogeneous systems 

[HiiniiHiiaiiiiiisiiniiiiiiisiEQiEiiEaEiiEiiEa 

for example, the relaxation of correlation functions in space and in time is generally 
in an exponential form, which defines a quench-dependent correlation length and a 
relaxation (or decoherence) time. Many basic features of the relaxation process can 
be successfully explained by a quasiparticle picture [33l [Ml |34] : after a global quench 
quasiparticle are created homogeneously in the sample and move ballistically with 
momentum dependent velocities. The behavior of observables in the stationary state 
is generally different in integrable and in non-integrable systems. For non-integrable 
models, thermalization is expected [H [H [H [H [H [IZl [H [H EOl El E2] and the 
distribution of an observable is given by a thermal Gibbs ensemble; however, in some 
specific examples this issue has turned out to be more complex [23l EH EH EI] ■ By 
contrast, it was conjectured that stationary state averages for integrable models are 
described by a generalized Gibbs ensemble [12| . in which each integral of motion is 
separately associated with an effective temperature. 

Concerning quantum quenches in inhomogeneous systems, there have been only a 
few studies in specific cases; for example, entanglement entropy dynamics in random 
quantum chains [35l ESI E2 and in models of many-body localization [38l ESj . In some 
of these cases the eigenstates are localized, which prevents the system from reaching 
a thermal stationary state. 

A special type of inhomogeneity, interpolating between homogeneous and 
disordered systems, is a quasicrystal [401141] or an aperiodic tiling [42]. Quasicrystals 
are known to have anomalous transport properties [43l[44], which is due to the fact that 
in these systems the long-time motion of electrons is not ballistic, but an anomalous 
diffusion described by a power law. One may expect that the quasiparticles created 
during the quench have similar dynamical behavior, which in turn affects the relaxation 
properties of quasicrystals. 

Quasicrystals of ultracold atomic gases have been experimentally realized 
in optical lattices by superimposing two periodic optical waves with different 
incommensurate wavelengths. An optical lattice produced in this way realizes a 
Harper's quasiperiodic potential [171[35], for which the eigenstates arc known to be 
either extended or localized depending on the strength of the potential. Different 
phases of Bose-Hubbard model with such a potential have been experimentally 
investigated [45l [46]. There have been also theoretical studies concerning the 
relaxation process in the Harper potential [49l [50] . 

In this paper we consider the nonequilibrium quench dynamics of the quantum 
Ising chain in one-dimensional quasicrystals. The quantum Ising chain in its 
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homogeneous version is perhaps the most studied model for nonequihbrium relaxation 
[nilSllMllMllMllillSTlliaililSlllHllIIllSllSllMl Our study 

extends previous investigations in several respects and seek to obtain new insights 
into quench dynamics in inhomogencous systems. We focus on the Fibonacci lattice, 
for which many equilibrium properties of the quantum Ising model are known 
[HlllMllMllIQllIIlllSlIZS]; to our knowledge this is the first study of quantum quenches 
in such a lattice. Using free-fermionic techniques |74j . we numerically calculate the 
time-dependence of the entanglement entropy as well as the relaxation of the local 
magnetization for large lattices. The numerical results are interpreted by a modified 
quasiparticle picture, in which the quasiparticles are represented by wave packets; we 
also obtain diffusive properties of the wave packets. 

The structure of the paper is as follows. The quasipcriodic quantum Ising model 
and its equilibrium properties are described in section [51 The global quench process 
and some known results for homogeneous and random chains are presented in section[3] 
Our numerical results for the quasipcriodic chain are presented and interpreted in 
section |4l This paper is closed by a discussion; some details of the free-fermionic 
calculation of the local magnetization arc presented in the appendix. 



2. The Model and its equilibrium properties 



We consider the quantum (or transverse) Ising model defined by the Hamiltonian: 



-4 



(1) 



where erf and af are Pauli matrices at site i. 
dependent, which are parameterized as: 



The interactions, J^, are generally site 

(2) 

where r > is the amplitude of the inhomogeneity, and the integers fi are taken from 
a quasipcriodic sequence. 

Quasipcriodic lattices can be generated in different ways, such as by the cut-and- 
projcct method. Here we use the following algebraic definition for a one-dimensional 
quasipcriodic sequence: 

1" 



/. = 1 



(3) 



where [x] denotes the integer part of x, and tj > 1 is an irrational number. The 
Fibonacci sequence generated by the substitution rule: — > 01 and 1 — > starting 
with corresponds to the formula in with the golden mean uj = (\/5 + l)/2. The 
parameter J in is fixed with J = r^P, where 



lim 



1 



1 



(4) 

1 represents the 



L^oo L OJ 

is the fraction of units 1 in the infinite sequence. Note that r 
homogeneous lattice. 

The model in ([!]) can be transformed into a free fermionic problem [74] (see 
Appendix). The spectrum of free-fermionic excitations, in (jA.ip . plays a key role in 
equilibrium and non-equilibrium properties of the system. In equilibrium and in the 
thermodynamic limit the model has a quantum critical point at h ~ he, the properties 
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of which arc controlled by the low-energy excitations. The value of is determined 
by the equation [75] In he = In J, where the overbar denotes an average over all 
sites. With the parameterization given above, the critical point is given hy he = 1, 
independently of r. The lowest gap, AE = ei, is zero for h < he, and it vanishes as 
AE ^ {he — h)" , as h approaches he. The singularity of the gap, measured by the 
gap-exponent v ~ 1, does not depend on r; the same is true for the singularity of the 
specific heat: Cy ^ ln\h — he\. Thus the transition belongs to the Onsager universality 
class irrespectively of r. This means that the quasiperiodic modulation of the couplings 
represents an irrelevant perturbation at the critical point of the homogeneous model 
[76] . For h < he the system is in the ordered phase, so that the local magnetization 
at site I is mi > 0. The bulk (surface) magnetization vanishes at the critical point as 
mb{h) ~ {he — hY^^ {mi{h) ~ {he — h)^/'^). For h > he the system is in the disordered 
phase and the local magnetization vanishes in the thermodynamic limit. 

While in equilibrium only the low-energy excitations are of importance, the 
complete energy spectrum contributes to nonequilibrium properties, which are 
investigated in this paper. 



3. Nonequilibrium properties of homogeneous and random chains 



We consider a quench process in which at time t = the strength of the transverse 
field is changed suddenly from ho to another value, say h. The initial Hamiltonian 

with ho for t < is denoted by "Ho, and its ground state is ^If^'^- For t > 
the new Hamiltonian H with h governs the coherent time-evolution of the system; 
for example an observable, represented by the operator A, has the time-evolution 
in the Heisenberg picture as: A{t) = exp{it'H)Aexp{—it'H), and its expectation 

value for i > is given by A{t) = (^i'lf^ A{t) ^'o'^''^ Dynamics of the system 
out-of-equilibrium is governed by the complete spectrum of H and not only by the 
lowest excitations. Therefore, Hamiltonians with different spectral properties will have 
completely different nonequilibrium properties. 

The form of the inhomogeneity in the couplings is generally crucial to the 
spectrum of a Hamiltonian. For example the spectrum of the homogeneous quantum 
Ising chain is absolutely continuous, thus all the eigenstates are extended. By contrast, 
the random chain has a singular point spectrum and the eigenstates are localized. The 
spectrum of quasiperiodic chains lies between the above mentioned two limiting cases 
[771175] ; for example, the spectrum of the Fibonacci chain defined in ^ is given by a 
Cantor set of zero Lebesgue measure, thus purely singular continuous [75] . 

Below we first briefly review nonequilibrium properties of the entanglement 
entropy and local magnetization after a quench in the homogeneous chain and in 
random chains. 



3.1. Entanglement entropy 

The entanglement entropy, Se{t), of a block of the first i sites in the chain is 
defined as Si{t) = Ti i[p i{t) In pi{t)] in terms of the reduced density matrix: pe{t) = 
Trj>£| !Z'o(i))( !^o(i)|- Here | 'I^o{t)) denotes the ground state of the complete system at 
time t > 0. Details of the calculation of Sg {t) in the free-fermion representation can 
be found in the appendix of [50] . 
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For the homogeneous chain (corresponding to the case with r = 1 in ([2])) in the 
limit L ^ CO and ioi l':$> 1 the results can be summarized as follows [331 131] '■ 

where Wmax is a maximum velocity. For a quench to a quantum critical point, the result 
in ([5]) is a consequence of conformal invariance [33]; for other cases, this behavior can 
be explained in the frame of a semiclassical (SC) theory [33l[34l: entanglement between 
the subsystem and its environment arises when two quantum entangled quasiparticles, 
which are emitted at t = and move ballistically with opposite velocities, arrive in 
the subsystem and in the environment simultaneously. The prcfactors a = a{ho,h) 
and /3 = P{ho,h) have been exactly calculated [81] and these agree with the results 
obtained from the SC theory [34]. In [82] a{ho = 0, h) has been evaluated in a closed 
formula, which is a continuous function of h but at the critical point /i = 1, its second 
derivative is logarithmically divergent. 

In the random chain the excitations are localized and therefore the dynamical 
entanglement entropy approaches a finite limiting value. When the quench is 
performed to the random quantum critical point, the average entropy increases ultra- 
slowly as log log t |36| . This behavior can be explained in terms of the strong disorder 
renormalization group [HS] [HU [36l [39] . 



3.2. Local magnetization 

Another quantity we consider is the local magnetization, mi(t), at a position, /, of 
an open chain. Following Yang [85] this is defined for large L as the off-diagonal 
matrix-element : 

mz(i) = (-Z'rKWI ^['^) , (6) 

where is the first excited state of the initial Hamiltonian. Calculation of the 

magnetization in terms of free fermions is outlined in Appendix. 

For the homogeneous chain the time-dependence of the local magnetization has 
been numerically calculated in [M] . For the quench performed within the ordered 
phases, Hq < 1 and h < 1, the results in the limit L ^ oo and I ^ 1 are given by: 

^ \ exp(-VC), t-»llv^^^ , ^ ' 

where the relaxation (decoherence) time r and the correlation length ^ depend on the 
quench parameters /iq and h. Exact expressions of these quantities have been derived 
recently [SU [Ml IM]- In the small /iq and h limit, accurate results can also been 
obtained from the SC theory [M] [55] . In this framework the quasiparticles in terms of 
the cr operators are represented by ballistically moving kinks. Each time when a kink 
passes site I, the erf operator changes sign; thus kinks that pass a site an even number 
of times have no effect on the local magnetization. Summing up contributions of all 
kinks, we obtain the functional form in ([7]). If the quench is performed close to the 
critical point, the kinks have a finite width; this effect can be taken into account in a 
modified SC theory [53] [Hi], which provides exact results. 

For quenches involving the disordered phase with /iq > 1 and /or ft, > 1 , the results 
obtained numerically [34l El] or analytically by the form-factor approach [61] [63] [64] 
indicate that the local magnetization changes sign during the relaxation process. The 
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period of these oscillations increases and becomes divergent as h ^ 1+, signaling a 
dynamical phase-transition in the system. 

In an disordered chain away from the random quantum critical point the bulk 
magnetization approaches a finite limiting value, which reflects the localized nature of 
the excitations. After a quench performed to the critical point, the average bulk 
magnetization has been found to vanish asymptotically in a very slow way |87j . 
mb{t) ^ \n~^{t), where A > is a disorder dependent constant. 



4. Results for quasiperiodic chains 

In this section we present our results for the quasiperiodic quantum Ising chain 
after a global quench, obtained by numerical calculations based on the free-fermion 
representation of the model. We concentrate on the Fibonacci chain with the 
parameter to defined in ([3]) being the golden mean. We consider finite chains with a 
length fixed to a Fibonacci number Fn . We have calculated the entanglement entropy 
and the local magnetization for system sizes up to L = Fu — 1597. Below we present 
results for these two quantities separately. 




Figure 1. Dynamical entropy after a quenchi from /iq = to various values of 
h at the aperiodicity parameters (a) r = 0.75, (b) r = 0.5 and (c) r = 0.25. 
The solid lines are the results for L = Fig = 987, and the dashed lines (only at 
h = 0.25, h = 0.5 and h = 0.75) correspond to the data for L = F17 = 1597. 

For a chain of total length Fn with periodic boundary conditions, we have 
calculated the entanglement entropy Sg between a block of length £ = Fn-2 and 
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its environment which has a length of -Fn-i- Various values of < r < 1 for the 
inhomogeneity amplitude were considered. We start our numerical calculations from 
the fully ordered state with /iq = to a state with h > both in the ordered and in 
the disordered phases, as well as at the critical point. The numerical results for 
Se{t) — 5^(0) are shown in figure [1] For all cases considered, Se{t) exhibits two 
time-regimes: in the late-time regime, the entropy is saturated to an L dependent 
value, similar to the behavior for the homogeneous chain; in the early-time regime, it 
increases with time as a power-law form: 

S{t) ^ t'' , (8) 

with some exponent cr < 1. The exponent <t depends on the value of the transverse field 
in the final state, while it does not vary (significantly) with the initial h^. The values 
of (J for r = 0.25, 0.5 and 0.75 are plotted in figure[5] for all cases considered, a reaches 
its maximum at the critical point /i = 1, and the increase with h in the ordered phase 
{h < 1) is much faster than the decrease in the disordered phase {h > 1). Furthermore, 
wc have found that the exponent a decreases with stronger inhomogeneity, that is with 
smaller value of r. 

The power-law time-dependence of the entanglement entropy in (|S]) is a new 
feature of the quasiperiodic system: the increase in entropy is slower than in the 
homogeneous chain, but faster than in a random chain. This behavior can be 
explained in terms of quasiparticles that are emitted at time t = 0, and subsequently 
move classically by anomalous diffusion which has a power-law relationship between 
displacement and time, x ^ . with a diffusion exponent < 13 < 1. The dynamics 
of the quasiparticles will be studied in more detail in section [4.31 



^.2. Local magnetization 

The local magnetization, mi{t), is calculated for open chains of length L = Fn- 
Generally mi{t) has a monotonous position dependence: mi^{t) > mi^{t) for li < 
I2 < L/2. We measured the magnetization at site I = Fn~i, which is considered as 
the bulk magnetization and denoted by mb{t). We have also studied the behavior of 
the surface magnetization, mi{t), for which some exact results are obtained. 

Wc study the asymptotic behavior of the surface magnetization (given in (jA.17|) ) 
for large t after a quench. If the quench is performed to the ordered phase, h < 1, 
the lowest excitation energy is ei = (i.e. cos(eit) = 1); consequently Pi,2fe-i(t) in 
(jA.Sp has a time independent part. This results in non-oscillating contribution to the 
surface magnetization by 

L 

m = W)Y.^^ij)^i^ij)^ (9) 

which corresponds to its asymptotic value. Recall that 3'i(l) = mi{h,t = 0) is the 
equilibrium surface magnetization |881 189| , which is finite for h < 1, and zero in the 

disordered phase. Similarly, $^°^(1) > for ho < 1 and zero otherwise. From this 
follows that the nonequilibrium surface magnetization is mi > 0, if both h < 1 and 
ha < 1. Otherwise the stationary surface magnetization vanishes. If the quench starts 
from the fully ordered initial state ho = 0, then 'I>i(j) = (5ij and rfTT = $i(l); thus 
we obtain the simple relation: 

mr(/i) = [mi(/i,t = 0)]2 , (10) 
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Figure 2. Double logarithm of tiic hulk magnetization as a function of tlie 
logarithm of the time. During the quench the transverse field is changed from 
ho = 0. to different values of h at the aperiodicity parameter r = 0.75 (panel 
(a)), r = 0.5 (panel (b)), r = 0.25 (panel (c)). The length of the chain is 
L = Fit = 1597 and the magnetization is considered at site I = -Fie = 987. 



which is generaUy vahd between the stationary value of the nonequihbrium surface 
magnetization and its equihbrium value. From (jlOp follows that the critical exponents 
are related as: /3"° = 2/3s. For the Fibonacci chain there is a linear variation close to 
the critical point: rrTi{h) ^ 1 ~ h^, thus = 1. 

We numerically calculated the time-dependence of the bulk magnetization after 
a quench from the fully ordered initial state, ho = 0, to different values of h. For fixed 
values of the inhomogencity r = 0.25, 0.5, 0.75, the results for the double logarithm of 
|mb(t)| are shown in figure [5] as functions of Int. In each case one can observe a linear 
dependence, which implies that the magnetization has asymptotically a stretched 
exponential time dependence: 

mb{t) - A{t)cxp{~Ctf') . (11) 

Here the prefactor A(t) changes its functional form at some critical value, h*{r), of the 
transverse field. For small transverse fields h < h*{r), the magnetization stays positive 
with A{t) > during the process. For larger transverse fields h > h*{r), A{t) changes 
sign in an approximately periodic fashion; thus in this region the magnetization 
oscillates in time; the oscillations are characterized by a typical time period, denoted 
by t*{h,r). As some critical value, /i*(r)+, of the transverse field is approached, 
h — ?> /i*(r)+, t*{h,r) becomes divergent. Thus at h*{r) there is a dynamical phase 
transition, similar to that noticed in the homogeneous model at ft-c = 1. Interestingly, 
for the quasiperiodic model this dynamical phase transition does not coincide with 
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the equilibrium phase transition. Estimates of h*{r) versus r are shown in figure [3] 
the data are well approximated by a power-law h*{r) ^ r" with a = 0.24(3) |90) . 




0.1 1 

r 



Figure 3. Position of tlie dynamical critical point for different values of the 
apcriodicity parameter in a double-logarithmic plot. The straight line has a slope 
a = 0.24. 

The exponent ^ describing the decay of the local magnetization dependents both 
on h and r; by contrast, it does not vary significantly with ho, at least for ho < h. Our 
results for the critical exponents fi = n{h,r) are plotted in figure [5] for r = 0.75, 0.5 
and 0.25 as functions of h. The exponent fi reaches its maximum at the dynamical 
phase transition point h* (r) . 

4-3. Interpretation by wave packet dynamics 

As known from previous studies on the homogeneous chain, dynamical features of 
the entanglement entropy and the local magnetization can be well described by 
the dynamics of quasiparticles. To understand the dynamical properties of the 
quasiparticles emitted after a quantum quench in the quasiperiodic lattice, we regard 
the quasiparticles as wave packets and study their dynamics using a method that has 
been applied for studies of transport properties of quasicrystals |441 [9Tj . 

First we notice that the time dependencies of the fermion operators for the 
quantum Ising chain are given in (jA.7[) and (jA.8|) . To construct a wave packet, which 
connects sites k and / at time t, we combine the four time-dependent factors in (jA.Sp 
in the form: 

WlMt) = I [P2l~1.2k-1 + P2l,2k + l{P2l-l,2k - P2L2k~l)] , (12) 

which is localized at t = since Wi,k{0) = di,k- The width of the wave packet starting 
from site k after time t is given by: 

i_ 

2 

(13) 



d{k,t) = 



j2ik-ir\w,Mt)\' 



. I 



One expects that for large t the wave packet in the infinite system shows anomalous 
diffusion in the form d{k,t) ~ t^^*^) with a diffusion exponent D{k), which generally 




depends on the starting position. After a global quench, quasiparticles are emitted 
everywhere in lattices, therefore d{k, t) should be averaged over different initial 
positions, 

d{t) = d(M) ■ (14) 

In our numerical calculations chains of length L = i^iy = 1597 with periodic 
boundary conditions were considered. First we have confirmed that the wave packet 
constructed in our method moves ballistically in the homogeneous chain (with r = 1), 
corresponding to Z? = 1. In the quasiperiodic chains the motion is indeed anomalous 
diffusive with Z3 < 1 , which is seen in figure |4] where the average widths of the wave 
packet are presented as functions of time in a log-log plot for various values of h and 
r = 0.75, 0.5 and 0.25. The diffusion exponent D for given h and r corresponds to 
the slope of the linear part of the function. 

The variation of D with ft, at a fixed r is shown in figure [SJ compared with 
the exponent a for the entanglement entropy and the exponent /i for the local 
magnetization. Here one can observe that the agreement between these three 
exponents is very good for h < h*{r), i.e. in the non-oscillating phase, but the 
exponent for the magnetization deviates in the oscillating phase (h > h*{r)). The 
discrepancy in the oscillating phase implies that the semiclassical picture breaks down 
in the oscillating phase, where the quasiparticles cannot well described by the moving 
kinks in the magnetization. 
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Figure 5. Scaling exponents calculated from the time-dependence of the width 
of the wave-packet, from the entanglement entropy and from the magnetization 
at different values of h for r = 0.75 (panel (a)), r = 0.5 (panel (b)), r = 0.25 
(panel (c)). The full lines connecting the diffusion exponents are guides for the 
eye. 



5. Discussion 

In this paper wc have studied the nonequihbrium dynamics of quasiperiodic quantum 
Ising chains after a global quench. In a quench process, the complete spectrum of 
the Hamiltonian is relevant for the the time evolution of various observables. For the 
quasiperiodic quantum Ising chain the spectrum is in a very special form, which is 
given by a Cantor set of zero Lebesgue measure, thus it is purely singular continuous. 
We have calculated numerically two quantities: the dynamical entanglement entropy 
and the relaxation of the local magnetization. The entanglement entropy is found 
to increase in time as a power-law (see ([HI)), whereas the bulk magnetization decays 
in a stretched exponential way (see (jlip V Both behaviors can be explained in a 
quasiparticle picture, in which the quasiparticles move by a anomalous diffusion in the 
quasiperiodic lattice. The diffusion exponent has been calculated by a wave packet 
approach, and good agreement has been found with the exponents that we obtained for 
the entropy and for the magnetization. We note that the anomalous dynamics found 
in the global quench process is similar to the transport properties of quasicrystals. 

Relaxation of the bulk magnetization is found to present a nonequilibrium 
dynamical phase transition. The non-oscillating phase, in which the magnetization 
is always positive, and the oscillating phase, in which the sign of the magnetization 
varies periodically in time, is separated by a dynamical phase transition point, at 
which the time-scale of oscillations diverges. This singularity point, due to collective 
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dynamical effects, is different from tlie equilibrium critical point. 

Similar nonequilibrium dynamical behavior is expected to hold for other 
quasiperiodic or aperiodic quantum models as long as the spectrum of the Hamiltonian 
is also purely singular continuous; there is a large class of such models, for example 
the Thue-Morse quantum Ising chain. If, however the spectrum of the Hamiltonian 
of the model is in a different type, such as the Harper potential which has extended 
or localized states, the nonequilibrium dynamics is expected to be different than the 
case we consider in this paper. 

Acknowledgments 

FI is grateful to D. Karevski, H. Rieger and A. Siito for discussions. FI acknowledges 
support from the Hungarian National Research Fund under grant No OTKA K75324 
and K77629; he also acknowledges travel support from the National Science Council 
(NSC), Taiwan, under Grant No. 101-2912-1-004-518. YCL is supported by the 
NSC under Grants No. NSC 98-2112-M-004-002-MY3 and No. NSC 101-2112-M- 
004-005-MY3; she also gratefully acknowledges support from the National Center for 
Theoretical Sciences of Taiwan. 

Appendix A. Pree-fermionic calculation of the time-dependent local 
magnetization 

Using the Jordan- Wigner transformation the Hamiltonian in ^ can be written in 
a quadratic form in fermion operators, which is diagonalized through a canonical 
transformation, yielding 



Here the energies of free fermionic modes, e^, are given by the solutions of the 
equations: 



with normalized vectors and ^k- For the calculation of the local magnetization we 
use free boundary conditions with Jl = 0, corresponding to ^k{L -I- 1) = Vl'fc(O) = 0. 

To calculate of the local magnetization in ([5]), we introduce at each site two 
Majorana fermion operators, 0,21-1 and 0,21, defined as 




fe=i 



(A.l) 



h<i>k(l)-Jk'^k{l + l) , 
Jk-i-^k{l-l)-h'^k{l) 



L 



a2i-i 



k=l 



L 




(A2) 



These satisfy the commutation relations: 

a-l^ = au {auhk} = 2(5;^^ . 
The spin-operators are expressed in terms of the Majorana operators as: 



(A3) 



.X 



(A.4) 



Nonequilibrium quench dynamics in quantum quasicrystals 



13 



and the local magnetization in (jB]) is then given as the expectation value of product 
of fcrmion operators with respect to the ground state. 

2/-1 

where we have used: 1 1?}*^-*) = 771 1 l/'g*'''). The expression in (jA.5[) - according to Wick's 
theorem - can be expressed as a sum of products of two-operator expectation values. 
This can be written in a compact form of a Pfaffian, which in turn can be evaluated 
as the square root of the determinant of an antisymmetric matrix: 



i~i 



■ {di{t)d2i-i{ti)) {di{t)r/i) 
■■ (a2(i)«2i-i(i)) {a2[t)m) 

(a2i-2(<)«2;-i(i)) {d2i~2{t)rii) 
{d2i-i{t)rii) 



= ±[detC,j]^/^ 



(A.6) 



where Cij is the antisymmetric matrix Cy = —Cji, with the elements of the Pfaffian 
(jA.6|) above the diagonal. (Here and in the following we use the short-hand notation: 

(...)=(^ri--.i^r>-) 

The time evolution of the spin operators follows from the time dependence of 

the fermion operators: Inserting ril{t) = e***^*??]^ and r]k{t) = e~^*'^'"qk into (|A.2p one 
obtains 

2L 

dra{t) = ^ Pra,n{t)dn , (A.7) 



with 



^2;-i,2fe-i = ^cos(e5i)$,(/)$,(fc), 



^2;-i,2fe 
^2;,2fc-i 

P2l.2k 



-^sin(e,0$«(0*9(fc) 
<? 

^sin(e,t)$,(fc)vI/,(0 , 

q 

^cos(e,t)vl/,(0*,(fc) . 



The two-operator expectation values are given by: 

fci ,k2 

The equilibrium correlations in the initial state with a transverse field ho are: 

(a2m-ia2n-l) = {d2md.2n) ~ ^m,n, 
{d2m-ld2n) = —{d2md2n-l) = , 

where the static correlation matrix Gm]n is given by: 



(A.8) 
(A.9) 

(A.IO) 
(A.ll) 



Nonequilibrium quench dynamics in quantum quasicrystals 
Then (|A.9p can be written in the form: 

{t)dn{t)) = Sm,n + l^7n,n{t) , 
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(A.12) 



with 



m — 1,2^2 



.(0) 



- 2i-1.2m 



^fci,fc2-^2i-1.2fci m— 1.2fc2 — 1 

['^fc2!fcl-^2i-l,2fci-1^2rn,2fc: 



fcl.fc; 



- 2i,2m-l — 



^fei,/c2-^2(-1.2fci -P2m,2fe2-1 
^ X] [^i2!fcl-^2i,2fe2-P2,n-l,2fci-l 



fci ,fc2 



- 2l,2r. 



-l-r2m,2fc2 



fci ,fc2 

-.(0) 



— ^fci,fe2^2i,2fci-P2Tn,2fc2-l 

In (|A.6[) there are also the contractions: 

n™ = (<Z^ria™(i)'7i|>^r> 



(A.13) 



(A.14) 



where 



Thus finally the square of the local magnetization is given by the determinant: 



(A.15) 



mf(t) 






ri,2 


ri.3 


ri_2;- 


1 


Hi 


-ri,2 





r2.3 


r2,2i- 


1 


n2 


-ri,3 


~r2,3 





r3^2i- 


1 


Ha 


ri.2i-i 











^21- 


-Hi 






-n2,-i 








(A.16) 



As a special case, the surface magnetization is expressed as: 

L 
L 

i=i 



(A.17) 
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